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ABSTRACT 

Cosmological simulations indicate that cold dark matter (CDM) halos should be triaxial. Validating 
this theoretical prediction is, however, less than straightforward because the assembly of galaxies is 
expected to modify halo shapes and to render them more axisymmetric. We use a suite of iV-body 
simulations to investigate quantitatively the effect of the growth of a central disk galaxy on the shape 
of triaxial dark matter halos. In most circumstances, the halo responds to the presence of the disk 
by becoming more spherical. The net effect depends weakly on the timescale of the disk assembly 
but noticeably on the orientation of the disk relative to the halo principal axes and it is maximal 
when the disk symmetry axis is aligned with the major axis of the halo. The effect depends most 
sensitively on the overall gravitational importance of the disk. Our results indicate that exponential 
disks whose contribution peaks at less than ~ 50% of their circular velocity are unable to modify 
noticeably the shape of the gravitational potential of their surrounding halos. Many dwarf and low 
surface brightness galaxies are expected to be in this regime, and therefore their detailed kinematics 
could be used to probe halo triaxiality, one of the basic predictions of the CDM paradigm. We argue 
that the complex disk kinematics of the dwarf galaxy NGC 2976 might be the reflection of a triaxial 
halo. Such signatures of halo triaxiality should be common in galaxies where the luminous component 
is subdominant. 

Subject headings: cosmology: theory — dark matter — galaxies: halos — halos: shapes — halos: 
structure — methods: numerical 



1. INTRODUCTION 



Cold dark matter (CDM) halos are found to be triax- 
ial in cosmolog ical A^-body si mulatio ns fe.g..lFrenk et al.l 
1988; Dubin ski fc Carlberd [T99ll : IJing fc Sutd 120021 : 
Bett et al.ll2007l ). This finding has motivated a number of 



studies designed to constrain halo shapes using a variety 
of probes, incl uding the gravitationa.1 lensing of distant 
galax ies (e.g., iHoekstra et al.l [2004t iMandelbaum et aTl 
|2006( ) ; the m orphology of tidal streams in the Milky Way 
(MW ) (e.g., Ilbata et al.l WM . Hfeh^ [200l ILaw et al.l 
20091) : the kinem atics of polar ring galaxies (e.g. 



Sackett fc Sparkc 1990); the fiaring of galactic disks (e.g.. 
Oiling fc Merrifield.2000 '); and the X-ray emission from 



hot gas in galaxies and clusters (e.g.. iKolokotronis et al.l 
[200lt IBuote "ei~all[200l . 

Despite these efforts, a clear picture of either conflict 
or agreement with CDM predictions has yet to emerge. 
This is in part due to the inherent difficulty of the obser- 
vational task, but also because most theoretical predic- 
tions rely on simulations where the effects of the bary- 
onic component of galaxies are neglected. The assembly 
of a central galaxy can modify substantially the shape 
of its surrounding dark halo (see, e.g., [P ubinsk i 1994.; 
Kazantzidis et al.ll2004atlAbadi et aIll20ia[Tlssera et al.1 



20101) by reshaping the box orbits that sustain tri axial- 
ity (e.g.. lDebattista et al.l [20081 : iValluri et al.ll2010t ). Al- 
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though there is reasonable consensus on the qualitative 
effects of baryons on halo shapes, there has been little 
quantitative work aimed at gauging the response of re- 
alistic triaxial halo models to the changes in the various 
parameters that characterize the central galaxy, such as 
its mass, size, or the timescale and mode of its assembly. 

Here we explore these issues using a series of controlled 
numerical simulations where a triaxial halo is evolved un- 
der the infiuence of a central disk galaxy. Our numerical 
experi ment s are similar in nature to those of iDubinskil 
()1994| ) and lDebattista et al.l (|2008[ ). However, we extend 
this work in several respects. For example, we focus on 
the shape of the gravitational potential rather than that 
of isodensity contours. The potential is less sensitive to 
the infiuence of substructures, which can induce substan- 
tial but transient changes in the local density whilst af- 
fecting little the overall potential. The gravitational po- 
tential is also a quantity of more direct relevance and ap- 
plicability to many observational studies of halo shapes. 

Further, we explore systematically several aspects of 
the growth of the galaxy that may in principle affect the 
shape of its surrounding halo. We consider not only the 
gravitational importance of the disk, but also its orien- 
tation relative to the halo principal axes, as well as the 
timescale and mode of its assembly. Lastly, we use halo 
parameters in agreement with the results of cosmologi- 
cal A^-body simulations and galaxy parameters which are 
consistent with observed scaling laws and span the wide 
range in surface brightness of observed galaxy disks. 

2. SIMULATIONS AND METHODS 

We investigate changes in the shapes of triaxial dark 
matter (DM) halos induced by the grow th of central disk 
gala xies modeled as rigid potentials fsee lVillalobos et all 
l200i^ for details). We measure the shape of the gravita- 
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Fig. 1. — Top and bottom-left panels: 
respectively, to the halo before and after the addition of the disk (dotted line). 

but its scale length varies: = 22.3, 8, and 2.8 kpc. Galaxy masses are consistent with the baryonic TF relation for a rotation speed 
of ~ 120 kms~^. The disk models span roughly two decades in surface density (or 5 mag in surface brightness) from values typical of 
LSB galaxies to intermediate (ISB) to "normal" (HSB) spirals. Radii are given in units of ij^j. The halo circular velocity profile peaks at 
Vmax ~ 120 kms~^ at rmax ~ 48 kpc. Bottom-right: Circular velocity pro file (top t wo cu rves) and disk contribution to the cir cular velocity 
(bottom curve) for NGC 2976. The disk contri bution is taken from Simon et al.l II2003I ). Circular velocities are taken from ISTmon et aLl 
||2p03) (filled circles) and 'Sp ekkens fc Sellwoocll (j2007) (open squares). The shaded area indicates the radius, r = 2.2R^, where the disk 
gravitational importance parameter rj = Vd/Vdrc is measured. For NGC 2976, ~ 0.7 kpc and rj ~ 0.5. 



tional potential by approximating the isopotential sur- 
faces by ellipsoids and calculating the principal axes a , 
b, and c (e.g., .Sprin gcl et al. 20M IHavashi et all 120071: 



lAbadi et al where a, 6, and c denote the major, 



intermediate, and minor axis, respectively. 

The triaxial DM halo models ar e constructed via suc- 
cessive mergers fMoore et al.l (|2004[) . The initial system is 
a 75, 000-particle, isotropic, sphericall y-syrn metric model 
generated by sampling thelH crnquist' (jl990( ) distribution 
function (IKazantzidis et aLl i2004b.) . A nearly prolate 
remnant results from a head-on collision between two 
such Hernquist models; more triaxial configurations are 
built by merging this remnant with other Hernquist mod- 
els on parabolic orbits of various inclinations relative to 
the principal axes of the first remnant. Each final triax- 
ial halo consists of 300, 000 particles. After the mergers 
are complete, we evolve the remnants for several dynam- 
ical times in order to ensure that equilibrium is achieved. 
The final configurations have mass distributions th at can 
be w ell approximated by NFW profiles (Navarr o et al.l 
\19M) (in the regions relevant to the present study) and 



are thus consistent with the results of cosmological CDM 
simulations. 

We adopt two different halo models in order to grow 
central disks. The first (model A) is a ne arly prolate 
merger remnant, with triaxiality parameter (jFranx et al.l 
IT99I T ={0? - b'^)l{a? - c2) - 0.8 for most radii. The 
second (model B) is strongly triaxial; T ^ 0.4 in the in- 
ner regions, increasing to T ^ 0.6 in the outskirts. These 
shapes are f airly typical of thos e found in cosmological 
simulations ()Havashi et al.ll2007D . The spin parameter A 
of halos A and B are 0.024 and 0.04 0, respectively , whic h 
are typical of cosmological halos fMacc io et al.l 120081) . 
The circular velocity and axial ratio profiles of these ha- 
los are presented in Figures [T] and [U respectively. 

We scale the parameters of halos A and B to match the 
mass and concentration of CDM halos with maximum 
circular velocity V^ax ~ 120 kms~^. For a halo of virial 
mass Mv 



5 X 10^^M(7) and concentration 



^ 9.3 

(jNeto et al.l 120071 ). the circular velocity profile peaks at 

:nax ~ 48 kpC. 

Disk galaxies forming in our triaxial halos must have 
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Fig. 2. — Axial ratios of the halo gravitational potential measured in the plane of the disk before (solid lines) and after the growth of 
a central galaxy, as a function of radius. Arrows indicate the radius (r = 2.2R^) where the disk contribution to the circular velocity is 
maximum. The dot- dashed, dashed, and dotted lines show halo axial ratios after the growth of the LSB, ISB, and HSB disk, respectively. 
Panels on the left show the case where the disk plane is perpendicular to the minor (z) axis of the halo, while panels on the right correspond 
to the case where the disk plane is perpendicular to the major {x) axis of the halo. Results are presented for = 10 Gyr but are insensitive 

to Td. 



rotation speeds comparable to Vmax in order to sat- 
isfy simultaneously the normalization of the Tully-Fisher 
(TF ) relation and the galaxy stellar mas s function (see, 
e.g.. iCroton et all 120061: iGuo et al.ll2010D. According to 
the b aryonic TF relation of iNoordermeer fc VerheiienI 
((200l . the mass of the luminous component of galax- 
ies with rotation speeds of order 120 kms^^ is Md ~ 
1.3 X IQ^^Mq. We model these galaxies as exponen- 
tial disks with three different values of the radial scale 
length, Rd = 2.8, 8.0, and 22.3 kpc. (For the disk verti- 
cal structure we assume an isothermal sheet with scale 
height Zd = 0.2 Rd, consistent with observations of exter- 
nal galaxies.) These values span the range of observed 
scale lengths for galaxies of this mass, from the very ex- 
tended low surface brightness (LSB) disks to intermedi- 
ate (ISB) to concentrated high surface brightness (HSB) 
"normal" spirals. 

Within the disk radius the circular speed of the 
disk-|-halo system is roughly the same for all models, ir- 
respective of the surface brightness/density of the disk. 
This is in agreement with observations, which show that 
the zero point of the TF relation is roug hly independent 
of surface brightness (Zwaa n et al.lll995 ). 

In addition to surface brightness, we have also explored 
the response of the halo to varying the timescale of disk 
assembly. This is accomplished by allowing the disk mass 
to grow linearly over three different timescales, rd=0 Gyr, 
1 Gyr, and 10 Gyr. The third choice ensures that the 
response of the halo is adiabatic, whereas the first, al- 
though unrealistic, provides a useful check of the sensi- 



tivity of our results to Td- We have also varied the mode 
of disk assembly by simulating cases where the mass of 
the disk grows in three discrete steps over 1 Gyr, in- 
stead of linearly. None of the results we obtain depend 
on either the choices for Td or the mode o f disk assem- 
bly (see also iBerentzen fc Shlosmanl I2006D , so we shall 
mainly concentrate on the "adiabatic" (t^ = 10 Gyr) ex- 
periments where the mass of the disk grows linearly with 
time. 

Because the precise alignment between the disk angu- 
lar momentum axis and the h ost halo principal axes is 
still a matter of debate (e.g.. iFaltenbacher et al.l 120051 : 
iZentner et al.ll2005l: iBailin et al.l l2007r the final param- 
eter we investigated is the orientation of the disk plane 
relative to the principal axes of the triaxial halo. For 
each halo, we placed the disk plane perpendicular to its 
major as well as its minor axis. In total, we performed 
24 simulations of the growth of a central disk galaxy in 
each of the triaxial halo models A and B, for a total of 
48 experiments. After the disk growth is complete, all 
simulations are continued for at least a further 1.5 Gyr 
in order to allow for equilibrium to be reached. In the ex- 
periments where the disk grows in three distinct events, 
one-third of the disk mass increases instantaneously at 
0.1, 0.6, and 1.1 Gyr. We only consider the final equi- 
librium configurations of the halos when measuring the 
shape of the potential. 

All numerical simulations were c arried out w ith the 
parallel TV-body code PKDGRAV (|Stadell [200l . The 
softening was set to e = 0.5 kpc. Each disk was modeled 
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Fig. 3. — Axial ratios of the inner (r < 30 kpc) potential, as a function of the disk gravitational importance parameter, rj = Vd/Vdicy 
measured at r = 2.2R^. Symbols and error bars indicate the average and the rms shape of the halo isopotential contours within 30 kpc, 
respectively. The leftmost point in each panel (r) = 0) corresponds to the halo before adding the disk, while the rightmost point {rj = 1) 
corresponds to a fully dominant axisymmetric disk. Intermediate values of rj show results for disks of different surface density: LSB, ISB, 
and HSB, from left to right. Results correspond to a disk growth timescale of = 10 Gyr. The orientation of the disk plane relative to 
the principal axes of the halo in each panel is indicated in the labels. For given rj, the effect is maximal when the disk symmetry axis is 
aligned with the halo major axis. The shaded area highlights the value of rj of NGC 2976. 



with N = 10^ static particles in order to ensure a smooth 
representation of the disk potential 

3. RESULTS 

Figure [2] shows the axial ratios of the halo isopotential 
contours measured in the plane of the disk, for various 
experiments. The growth of the disk modifies the halo 
shape, making it more axisymmetric. Halo A is nearly 
prolate so placing the disk plane perpendicular to the 
major axis has little effect on the shape of the 2D poten- 
tial on the disk plane. When the disk plane contains the 
major axis the potential is strongly non- axisymmetric, 
but becomes rounder after the disk is added. The effect, 
however, is minor for the case of the LSB and ISB disk; 
only the HSB galaxy is able to modify substantially the 
halo shape, increasing the inner axial ratio from ~ 0.6 to 
- 0.8. 

Before adding the disk, halo B is rather triaxial, with 
h/a ~ c/h ~ 0.8 in the inner regions. In the case of the 
HSB disk, placing the disk plane perpendicular to the 
halo minor axis renders the halo potential almost per- 
fectly axisymmetric (measured in the plane of the disk). 
On the other hand, the LSB and the ISB galaxies are 



again barely able to modify the shape of the halo. When 
the disk plane is perpendicular to the halo major axis 
the results are similar, although in this case the halo 
response to the ISB and HSB galaxies are comparable. 
Except very near the center, the potential remains far 
from axisymmetric in all cases. We note that the change 
in shape in the case of the HSB disk is noticeable out 
to ~ 30-50 kpc, well outside the region where the disk is 
gravitationally important. 

Figure |3] shows how the 3D halo shape changes in re- 
sponse to the disk growth. The axial ratios of the in- 
ner gravitational potential are presented as a function 
of the peak contribution of the disk to the circular ve- 
locity. This is measured by the parameter 77 = Vd/Vcitc, 
computed at r = 2.2Rd, the radius where the exponential 
disk contribution to the circular velocity is maximal. The 
shape of the potential is not exactly constant in the in- 
ner regions (see Fig. [2]) , and therefore we show the shape 
of the potential averaged inside 30 kpc. The trends we 
discuss are robust to reasonable changes in the averaging 
procedure. 

Results are presented for the "adiabatic" {tcl = 10 Gyr) 
disk growth and for both orientations of the disk plane 
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Fig. 4. — Axial ratios of the inner (r < 30 kpc) potential, computed in the plane of the disk, as a function of tj. Symbols and error bars 
are as in Figure \3\ Filled symbols show the shape of the halo isopotential contours; open symbols correspond to the total (disk+halo) 
potential. Circles and squares show cases where the disk growth timescale is T;j=0 Gyr and t^=10 Gyr, respectively. The orientation of 
the disk plane relative to the principal axes of the halo in each panel is as in Figure [2] 



relative to the halo principal axis, rj varies from ^0.25 
for the LSB galaxy to ~ 0.7 for the HSB model. Fig- 
ure |3] demonstrates that in most circumstances the disk 
growth renders the 3D potential of the halo more spher- 
ical. Because halo accelerations are minimal along its 
major axis, whereas the disk gravity is strongest along 
its symmetry axis, the effect of a disk on a triaxial halo 
is maximal when the disk symmetry axis is aligned with 
the halo major axis, and minimal in the opposite case. 

Figures [2] and |3] show that the modifications to halo 
shapes induced by the disks are relatively minor, except 
for the HSB galaxy. The response of the halo depends 
principally on the disk gravitational importance. These 
results seem insensitive to the timescale and mode of 
disk assembly. Figure |4] shows this explicitly. This figure 
shows the axial ratio of the inner potential as a func- 
tion of rj computed in the plane of the disk and confirms 
the conclusions advanced above. The change in shape is 
monotonic with rj and is relatively minor for the LSB and 
ISB models. Only the HSB galaxy is able to modify the 
halo potential noticeably. The transition seems to occur 
roughly at 77 w 0.5. Disks that contribute less than about 
50% of its circular velocity should be surrounded by halos 
that have preserved their original triaxiality and, there- 
fore, should in general show signatures of departures from 



axisymmetry. This is the regime of many LSB and dwarf 
galaxies, and it is therefore important to consider the 
possibility that their halos may be triaxial when an alyz- 
ing their kinematics. IBerentzen fc Shlosmanl (j2006D also 
concluded that maximal disks probably reside in nearly 
axisymmetric halos. 

4. DISCUSSION 

The sphericalization of DM halos by central galaxies 
has important implications for the interpretation of ob- 
servational constraints. For example, a nearly spherical 
halo is favored by the weak precession of the Sagittar- 
ius stream, a result that would be difficult to understand 
in the CDM paradigm unless the MW hal o is an out- 
lier i n the distribution of h alo shapes (e.s... Ilbata et alj 
120011: iMaiewski et~all I2003D . Indeed, given the gravi- 
tational importance of the MW disk (rj 0.8 for the 
canonical values of Vrot — 220 kms^^, Rd = ?> kpc, and 
Md = 5 X IO^^Mq), the MW halo should have been sig- 
nificantly affected by the disk. Nearly spherical halos 
may also a ffect the dynamics of bars and stars in galactic 
disks (e.g., Ildeta fc Hozumil 120001: lEl-Zant fc Shlosmanl 
[200l . and can inhibit the fueling of nuclear supermas- 
sive black holes by suppressing the large-scale asym- 
metries responsible for the loss of angular momentum 
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in the gas compone nt (e.g., iMerritt &: Quinlanl 119981 : 
ICallegari et aI1[20lol) . 

On the other hand, our results suggest that LSB 
halos should retain the triaxiality predicted by cos- 
mological simulations. This may help to explain the 
rathe r elongated shapes of the ultra-faint MW satel- 
lites (jMartin et al.l |2008() (although these can also be 
attributed to tidal effects in the gravitational field of 
the MW) and oddities in the kinematics of some dwar f 
spheroidal galaxies (see, e.g., iPenarrubia et al.l 120101 ). 
Many dwarf galaxies are also in the regime rj ^ 0.5 
where they are unable to modify their surrounding halos. 
Gaseous disks in such galaxies should exhibit departures 
from axisymmetry, unless the halos are either prolate or 
oblate and the disk plane coincides with that where the 
2D potential is axisymmetric. 

This could indeed be the case in nearly oblate ha- 
los, given the preference o f the angular mo mentum to 
align with the minor axis (jBett et al.ll2007| ). However, 
in nearly prolate halos (a more common occurrence ac- 
cording to A^-body simulations) disks whose angular mo- 
mentum aligns with the minor axis would feel a non- 
axisymmetric 2D gravitational potential. Non-circular 
motions should therefore be fairly common in the gaseous 
disks of dwarfs, and they could in principle be used to 
gauge the triaxiality of their surrounding halos. 

To first order, a gaseous disk in the non-axisymmetric 
potential of a triaxial halo would behave just like gas in 
a barred potential where the pattern speed of the bar is 
zero. A subdominant disk in a triaxial potential would 
thus exhibit the non-circular dynamical signature of a 
(slow) bar but with no obvious bar in the luminous dis- 
tribution. 

Interestingly, there is one system where all these con- 
ditions are met. NGC 2976 is a nearby dwarf spiral 
galaxy whose baryonic disk is subdominant, as shown 
in the bottom-right panel of Figure [TJ This figure shows 
that the contribution of the baryonic component peaks 
at about one-half of t he circular v elocity at r = 2.2Rd 
and therefore « 0.5. iSimon et al.l ([2 00 SI show that the 
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kinematics of the gaseous disk in NGC 2976 is highly 
complex, exhibiting large non-circular motions near th e 
center. These, according to lSpekkens fc SellwoodI (|2007| ). 
are best understood as the characteristic kinematic 
asymmetries imposed by a n m = 2 bar mod e in the gravi- 
tational potential (see also lHavashi fc N"avar ro 20061 . On 
the other hand, NG C 2976 has no obvious bar, at leas t in 
the optical (but see lMenendez-Delmestre et al.ll2007l ). so 
ascribing the origin of the non-circular motions to halo 
triaxiality is clearly tempting. 

If this interpretation is correct, then it would be sur- 
prising if other galaxies with subdominant baryonic com- 
ponents did not also show signs of being embedded in 
triaxial potentials. Indeed, one may even argue that the 
absence of such signatures in a significant fraction of un- 
barred LSB and dwarf galaxies would be quite difficult 
to accommodate within the standard CDM paradigm. 
Definitive conclusions on these issues require more so- 
phisticated theoretical modeling of the formation of sub- 
dominant disks in triaxial halos. Note that our models 
neglect, for example, the response of the disk to the triax- 
ial forcing of the halo, as well as a realistic accounting of 
the distribution of disk orientations relative to the prin- 
cipal axes of the halo. Nevertheless, our results suggest 
that a careful search for signatures of halo triaxiality in a 
statistically significant sample of dwarf and LSB galaxies 
would be warr anted. Steps in this d irection such as those 
taken by Tra chternach et all ()2009( ) should certainly be 
encouraged. 
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